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ABSTRACT 

We present SDSSJ092712. 65+294344.0 as the best candidate to date for a recoiling supermassive 
black hole (SMBH). SDSSJ0927+2943 shows an exceptional optical emission-line spectrum with two 
sets of emission lines: one set of very narrow emission lines, and a second set of broad Balmer and 
broad high-ionization forbidden lines which are blueshifted by 2650 km s~^ relative to the set of narrow 
emission lines. This observation is most naturally explained if the SMBH was ejected from the core 
of the galaxy, carrying with it the broad-line gas while leaving behind the bulk of the narrow-line 
gas. We show that the observed properties of SDSSJ0927-f 2943 are consistent with predictions and 
expectations from recent numerical relativity simulations which demonstrate that SMBHs can receive 
kicks up to several thousand kms~^ due to anisotropic emission of gravitational waves during the 
coalescence of a binary. Our detection of a strong candidate for a rapidly recoiling SMBH implies that 
kicks large enough to remove SMBHs completely from their host galaxies do occur, with important 
implications for models of black hole and galaxy assembly at the epoch of structure formation, and 
for recoil models. 

Subject headings: galaxies: active - galaxies: evolution - galaxies: individual 

(SDSSJ092712.65+294344.0) - quasars: emission lines 



1. INTRODUCTION 

The merging of two galaxies will produce a binary 
black hole at the center of the newly formed galaxy. If 
the two black holes do not stall, they will ultimately 
merge due to emission of gravitational wave radiation. 
The gravitational waves carry away linear momentum, 
causing the centre of mass of the coalescing BH system 
to recoil in the opposite direction (Peres 1962, Beken- 
stein 1973). Early analytical calculations predicted that 
mergers of non-spinning black holes can attain kicks with 
velocities of up to a few hundred kms~^ (e.g., Fitchett & 
Detweiler 1984, Favata et al. 2004, Blanchet et al. 2005, 
Damour & Gopakumar 2006), recently confirmed by nu- 
merical simulations (e.g.. Baker et al. 2006, Herrmann 
et al. 2007a, Gonzalez et al. 2007a). These velocities are 
above the escape velocity of dwarf galaxies, low-mass spi- 
rals, and high-redshift dark matter halos. If many BHs 
were removed from their hosts in the early history of 
the universe, this would have profound consequences for 
galaxy assembly and BH growth in the early universe, 
and would give rise to a population of interstellar and 
intergalactic BHs (e.g., Madau et al. 2004, Merritt et 
al. 2004, Madau & Quataert 2004, Haiman 2004, Yoo & 
Miralda-Escude 2004, Volonteri & Perna 2005, Volonteri 
& Rees 2006, Libeskind et al. 2006). 

Recent numerical relativity simulations of certain con- 
figurations of merging, spinning BHs have produced 
much higher recoil velocities, up to several thousand 
kms~^ (Campanelli et al. 2007a, b, Gonzalez et al. 
2007b, Tichy & Marronetti 2007, Herrmann et al. 2007b, 
Dain et al. 2008, Schnittman et al. 2008), scaling to an 
expected maximum around 4000 kms~^ (Campanelli et 
al. 20007a, b. Baker et al. 2008) for maximally spinning 
equal-mass binaries with anti-aligned spins in the orbital 
plane. These kick velocities exceed the escape velocities 
of even massive elliptical galaxies (Fig. 2 of Merritt et 
al. 2004) and therefore the new results reinforce and en- 



hance consequences studied earlier for the smaller kicks, 
with potentially far-reaching implications for the early 
phases of BH growth from early stellar-mass precursors 
or later intermediate-mass precursors (Schnittman 2007, 
Volonteri 2007) and consequently for the frequency of 
gravitational wave signals detectable with LISA (Sesana 
2007), for the scatter in the M — a relation (Libeskind et 
al. 2006), and for the offsets and oscillations of recoiling 
BHs in galaxy cores (Gualandris & Merritt 2008). 

The recoiling black holes will carry a fraction of nu- 
clear gas and stars with them (Merritt et al. 2004, 2006, 
Madau & Quataert 2004, Loeb 2007). They would be 
detectable spatially in the form of Seyfert or quasar ac- 
tivity offset from the galaxy core (Madau & Quataert 
2004), or in the form of broad emission lines kinemati- 
cally offset from the narrow emission lines (Bonning et 
al. 2007, Komossa et al. 2008). Because of the broad as- 
trophysical implications, the search for and actual iden- 
tification of such recoiling black holes is of great interest, 
and will place important constraints on BH growth dur- 
ing the epoch of structure formation, on predictions of 
maximum recoil velocity, and on arguments suggesting 
that the BH spin configurations leading to maximal recoil 
velocities should be rare in gas-rich mergers (Bogdanovic 
et al. 2007). 

Bonning et al. (2007) searched for recoiled SMBHs in 
the Sloan Digital Sky Survey (SDSS) database, looking 
for systematic kinematic offsets between broad-line gas 
attached to the recoiling BH, and narrow-line gas left 
behind. They did not find any promising candidate, and 
concluded that SMBH recoil with large kick velocities is 
relatively rare. 

Here, we present the best candidate to date for a 
recoiling SMBH, the quasar 

SDSSJ092712.65-f294344.0 (SDSSJ0927-H2943 here- 
after). Its unusual emission-line spectrum matches key 
predictions from the recoiled-SMBH scenario. We use 
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a cosmology with Hq^70 kins ^ 
r2A=0.7 throughout this Letter. 



Mpc" 



rjM=0.3 and 



TABLE 1 
Emission-line properties of SDSSJ0927+2943. 



2. OBSERVATIONS OF SDSSJ0927+2943 
2.1. Optical spectroscopy 

SDSSJ0927+2943 at redshift z=0.713 is a luminous 
quasar, observed in the course of the SDSS (Adelman- 
McCarthy et al. 2007), and was found by us in a system- 
atic search for active galactic nuclei (AGN) with high 
[OIII] velocity shifts. The SDSS spectrum, corrected for 
the Galactic reddening of E(B-V) — 0.021 mag, is dis- 
played in Fig. 1. The underlying continuum spectral 
energy distribution (SED) was modeled as a powerlaw 
with a best- fit slope of a ~ —0.1 (where f^y oc i^"). Each 
emission line was fit by a single Gaussian except the Fell 
multiplets, which were modeled by templates built from 
I Zw 1 (A > 3500 A, Veron-Cetty et al. 2004; A < 3500 
A, Tsuzuki et al. 2006). The redshifts of the Fell fines 
were tied either to Mgll (the UV multiplets) or to broad 
H/3 (the optical multiplets). 

Two systems of strong emission lines can be identi- 
fied in the spectrum, which we refer to as the "red" 
(r) and "blue" (b) systems. The red system con- 
sists of very narrow emission lines (red NELs, r-NELs 
hereafter) of [OIII]5007, [OII]3727, [NeIII]3869, faint 
[NeV]3426 and Balmer lines, all of them almost unre- 
solved (FWHM,obs([OIII]) = 230 kms^^; the widths of 
the narrow lines are all very similar, and were there- 
fore all fixed to the same value in order to derive 
fluxes). The blue system shows classical broad Balmer 
and MgII2798 emission lines (BELs), plus unusually 
broad NELs (blue NELs, b-NELs hereafter). AU hues 
of the blue system are blueshifted by about 2650 kms~^ 
relative to the r-NELs^ (see Tab. 1 for redshifts; the 
value of 2650 kms~^ is the shift between broad H/3 
and r-[OIII]). The b-NELs show broad [NeV] with a 
width of FWHM([NeV])==2080 kms^^, and broad [Nelll] 
with FWHM([NeIII])=1020 kms^i. [OIII] and [Oil] are 
present, too, with widths of 460 kms^^. The BELs ap- 
pear in Balmer fines and in Mgll with FWHM(H/3)=5740 
kms-i and FWHM(MgII)=3530 kms'^ Line ratios 
indicate AGN-like excitation in both systems, r-NELs 
and b-NELs. Emission-line properties are summarized 
in Tab. 1. 

In order to see whether any (faint) broad-line emission 
is also accompanying the r-NELs, we have performed the 
following test. We have first subtracted the best-fit con- 
tinuum. Fell multiplets, and NELs from the observed 
spectrum, and then fit the H/3 regime with two broad 
Gaussians. The redshift of the second Gaussian was fixed 
to that of the r-NELs, and its width constrained to be 
in the range 1000 — 5000 km s^^. No successful fit could 
be obtained if a contribution of the second Gaussian is 
enforced; its contribution is always consistent with zero. 
We therefore conclude that only the b-NELs are accom- 
panied by broad- line gas at the same redshift. 

The broadband SED of SDSSJ0927+2943 is rather 
blue with SDSS magnitudes of u ^ 18.71 ± 0.02, g = 
18.36 ± 0.02, r ^ 18.40 ± 0.02, i ^ 18.43 ± 0.02, 'z = 

^ Since the r-NELs arc all very narrow and arc all consistent 
with the same redshift, we assume henceforth that they define the 
restframe of the system. 



system 


line 


2 


FWHM= 


line-ratio 


r-NEL 


[OIII] 5007 


0.71279 


170 


10.1 




H/3 


0.71279^ 


170^ 


i.ot- 




H7 


0.71279^ 


170^ 


0.4 




NeIII]3869 


0.71279^ 


170^ 


0.7 




Oil] 3727 


0.71279^ 


170^ 


2.6 




NeV 3426 


0.71279^ 


170^ 


0.3 


b-NEL 


OIII 5007 


0.69713 


460 


6.7 




H/3 


0.69713^ 


460^ 


l.O'' 




NeIII]3869 


0.69678 


1020 


1.8 




Oil] 3727 


0.69801 


460^ 


1.5 




NeV] 3426 


0.69709 


2080 


4.0 


BEL 


H/3 


0.69770 


5740 


30.7 




H7 


0.69970 


3880 


10.7 




H5 


0.69996 


3260 


3.8 




Mgll 2798 


0.69832 


3530 


29.5 



^ fixed'' normalized to 
tion 



l.C^ corrected for instrumental resolu- 



18.31 ± 0.02, and GALEX^ (Martin et al. 2005) magni- 
tudes of NUV = 18.57 ± 0.06 and FUV = 19.49 ± 0.15. 
Assuming that the standard broad-line region (BLR) 
scaling relations (Kaspi et al. 2005) hold, we expect a 
BLR radius of ~0.1 pc and estimate a SIVIBH mass of 
SDSSJ0927-h2943 of Mbh = 6 10« IVIg from the width of 
up and the luminosity at 5 100 A. 

2.2. X-ray detection 

We have searched the X-ray archives for observations of 
SDSSJ0927-t-2943. The quasar is serendipitously located 
close to the edge of two ROSAT HRI images. The deeper 
exposure, of 19 ks duration, was performed in April-lMay 
1995. X-ray emission from SDSSJ0927+2943 is detected 
with a countrate of 0.0037±0.0007cts/s, which translates 
into a soft X-ray luminosity of Lx = 5 10^** erg/s (as- 
suming an X-ray spectrum with no intrinsic absorption, 
and with photon index Fx — —2.5). The second HRI 
image was taken in November 1994 with a duration of 
10 ks. SDSSJ0927+2943 is detected with a countrate of 
0.0047±0.0012 cts/s, consistent with the other measure- 
ment. The X-ray detection demonstrates the presence 
of an inner accretion disk (relevant for the discussion in 
Sect. 3.2). 

3. DISCUSSION 

3.1. The most unusual quasar pair known? 

Several pairs of quasars have been detected at pro- 
jected separations of 3-10". While some of them can 
be explained by lensing, others very likely represent real 
pairs (e.g., Kochanek et al. 1999; review by Komossa 
2003). 3 Is SDSSJ0927+2943 actually a binary quasar 

2 http://galcx.stsci.edu/GR2/ 

•^ In this context, it is interesting to mention HE0450-2958, a 
system consisting of an ultraluminous infrared galaxy (ULIRG) 
and a narrow-line Seyfert 1 (NLSl) galaxy. The host galaxy of the 
NLSl was not detected in HST imaging, suggesting it was a 'naked 
quasar' (Magain et al. 2005). Merritt et al. (2006) then revised 
downward the SMBH mass of the NLSl by a factor 10; the non- 
detection of the host galaxy is then consistent with the observed 
upper limit. Merritt et al. (2006) argued that the emission-line 
properties of the NLSl are incompatible with the scenario that the 
NLSl was actually ejected from the ULIRG, whether the mech- 
anism of ejection was gravitational wave recoil or three-body in- 
teraction of a binary BH with a third BH (Haehnelt et al. 2006, 
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Fig. 1.— SDSS spectrum of SDSSJ0927+2943, showing two sets 
of emission lines separated by a velocity of v ^ 2650 kms~^. Red: 
red set of narrow emission lines (r-NELs), blue and light blue: blue 
set of emission lines (b-NELs and BELs, respectively), grey: Fell 
spectrum. 

? We temporarily refer to this hypothetical system as 
SDSSJ0927+2943A,B. The difference in velocity of the 
two sets of emission lines exceeds the peculiar velocities 
observed in galaxy clusters, and is too large for the two 
galaxies to form a bound merger. Their redshift differ- 
ence corresponds to '^100 Mpc, if their redshifts are cos- 
mological. Consequently, we would then have a very un- 
likely projection effect, including not just one, but two in- 
trinsically extremely unusual AGN: SDSSJ0927-I-2943A 
with exceptionally broad neon lines of [Nelll] and [NeV] . 
And SDSSJ0927+2943B as one of the rare type-2 quasars 
with, in addition, exceptionally narrow emission lines 
[their observed width of ^230 kms^^ is below that typ- 
ically observed in quasar narrow-line regions (FWHM 
« 400 - 800 kms"^ Zakamska et al. 2003)]. This rare 
source would have to be projected by chance behind the 
other unusual AGN. 

We are therefore let to consider the alternative hypoth- 
esis that we actually see only one AGN, its SMBH and 
bound emission-line region having been ejected from the 
core with a line-of-sight speed of 2650 kms~^, while the 
bulk of the narrow-line region (NLR) and other interstel- 
lar medium (ISM) was left behind and shines in narrow 
emission lines. This is the scenario actually predicted 
in the context of the recent black hole recoil simulations 
discussed in Sect. 1, and it is the observational signature 
Bonning et al. (2007) searched for, and Komossa et al. 
(2008) paid attention to, but did not detect. 

3.2. A recoiling SMBH in SDSSJ0927+2943 

In this picture, SDSSJ0927-I-2943 underwent a merger 
in the past. Upon merging, its central SMBH recoiled 
taking with it the BLR and perhaps other very high ion- 
ization gas, leaving behind the NLR. Gas with a velocity 

Hoffman & Loeb 2006). Finally, Kim et al. (2007) pointed out 
that the ultraluminous IR emission is more likely associated with 
the NLSl galaxy. 



larger than the recoil velocity will remain bound to the 
SMBH, and a trail of partially bound gas may form. Ac- 
cretion activity may have switched off temporarily in the 
course of the binary SMBH merging, once the orbital de- 
cay time due to emission of gravitational waves had be- 
come smaller than the viscous timescale of the disk (e.g., 
Liu et al. 2003, Milosavljevic & Phinney 2005). After co- 
alescence, the inner disk will re-form quickly (Loeb 2007) 
and the accreting SMBH will then illuminate the bound 
gas, a trail of partially bound gas and swept- up ISM, the 
surrounding ISM/halo and the left-behind bulk of the 
NLR gas. 

The NLR gas left behind will retain memory of 
the original ionization from the pre-merger accretion 
phase for only limited time, given by the light-travel 
time, and the hydrogen recombination timescale, tj-cc ~ 
130 T^-^n'^^ yr, where T4 is the gas temperature in 10** K, 
and 713 is the gas density in units of 10^ cm~^. [OIII] will 
fade away more quickly than hydrogen lines (Binette & 
Robinson 1987). However, once the recoiling SMBH re- 
forms its inner accretion disk, it will re-illuminate parts 
of the NLR left behind (from larger distance than before) , 
and any surrounding ISM. Emission-line ratios would de- 
pend on the average density which would be higher in the 
classical NLR than in other ISM/halo gas, and on metal 
abundances which will be lower in the halo. 

Which aspects of this scenario do we observe in 
SDSSJ0927-h2943 ? The blue system of emission-Hnes, 
the BELs and b-NELs, represents gas bound to the re- 
coiling SMBH, seen in broad Balmer lines, broad Mgll 
and broad forbidden lines, especially [NeV]. After coa- 
lescence and recoil, matter orbiting the SMBH with a 
velocity much larger than the recoil velocity will remain 
bound to the SMBH (Merritt et al. 2006); i.e. matter 
within a region whose size is given by r^, < GMbh^"^- In 
our case, rb < 10^^ cm, which is about a factor of 3 larger 
than the BLR of SDSSJ0927-I-2943. The width of [Nelll] 
(of the b-NEL system) of 1020 kms~^ is too narrow for 
this gas to be originally bound to the SMBH (except 
in case of projection effects), but as the disk keeps ac- 
creting onto the recoiling SMBH, it will spread radially, 
and will also drive some outflows. The semi-broad [Oil] 
and [OIII] most likely have this same origin, and/or are 
from swept-up ISM in front of the path of the recoiling 
SMBH^. 

The red system of emission lines, the r-NELs, repre- 
sents the NLR gas left behind, and/or ISM surrounding 
the recoiling SMBH, not bound to it, but illuminated by 
its accretion disk. All r-NELs have the same profiles and 
are very narrow, and the emission-line ratios imply an 
AGN ionizing continuum. Can we distinguish between 
NLR and other ISM/halo gas, in terms of line ratios and 
line profiles ? Line ratios would depend on density and 
distance. As the recoiling SMBH moves away from the 
galaxy core, illumination of the NLR from an, on average 
larger, distance would decrease the degree of ionization, 
while illumination of low-density ISM surrounding the 
recoiling SMBH would increase the degree of ionization 
of the line-emitting gas. It is possible that we see a mix 

* Based on photoionization models (e.g., Komossa & Schulz 
1997) we expect little [Oil] emission from the gas bound to the 
recoiling SMBH; even though shielding geometries might be con- 
structed in which this actually becomes possible. 
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of these two processes. We do not expect to see pure 
very low-density halo gas. Its high ionization parameter 
would lead to a much higher degree of ionization than 
we actually observe, as we have also verified by test cal- 
culations with Ferland's photoionization code Cloudy (S. 
Komossa et al., in preparation). The density-sensitive 
[SII] doublet is outside the observed wave band, but will 
be accessible with NIR spectroscopy. 

Independent evidence for the origin of the r-NELs 
conies from their profiles which are very narrow. Their 
width is below the stellar velocity dispersion of cr,=260 
kms^^ predicted from the Mbh — cr* relation (Ferrarese 
& Ford 2005). This indicates that the r-NEL gas illu- 
minated by the recoiling SMBH either does not feel the 
full bulge potential, or the local velocity field of the illu- 
minated gas does not reflect the full velocity dispersion. 
This latter situation is expected if the recoiling SMBH 
was in a disk where the velocity pattern is predominantly 
rotational, so that we only see a fraction of the full pat- 
tern. 

The whole parameter space of BH recoil velocities is 
still being explored (e.g., Baker et al. 2006, 2007, Cam- 
panelli et al. 2007a, b, Choi et al. 2007, Gonzalez et 
al. 2007a, b, Herrmann et al. 2007a, b, Pollney et al. 
2007, Schnittman et al. 2008, and references therein; see 
Pretorius 2007 for a review). Several recent calculations 
focused on (almost) equal mass BHs with specific spin 
configurations such that kick velocities are maximized. 
The line-of-sight velocity we measure is comparable to 
the recoil velocities predicted by the runs of Gonzales et 
al. (2007b; Wkick « 2500-2650 km s^^) and by Dain et al. 
(2008; 3300 kms~^). The scaling formulae of Campanelli 
et al. (2007a) and of Baker et al. (2008) predict maximal 
recoil velocities of ^^3800 kms~^. In order to reach a high 
kick velocity, the pre-merger BHs of SDSSJ0927-I-2943 
must have been rapidly spinning and of nearly equal 
mass; i.e. the galaxy hosting SDSSJ0927-I-2943 must 
have undergone a major merger. 

Assuming that the recoiling SMBH carried with it an 
amount of mass which is not larger than its own mass, 
and that most of that gas is ultimately available for ac- 



cretion, we can estimate an upper limit on the duration 
of the quasar activity, iq. We further assume that the gas 
continues to accrete at its current rate of L/L^dd — 0.1, 
where L^dd is the Eddington luminosity, and we use a 
radiative efficiency of 77=0.37 which is appropriate for 
rapidly spinning BHs. This implies an upper limit on 
the lifetime of quasar activity of tq « 10^ yr. If the 
space velocity of the recoiling SMBH was close to the 
maximum possible velocity, it would reach a projected 
separation from the core on the order of a few kiloparsec 
within '^10^ yr. 

Future HST imaging of SDSSJ0927-f 2943 wiU be valu- 
able to distinguish whether the host galaxy shows a dis- 
torted morphology as would be expected if the merger 
was recent. With its spatial resolution of 0.1", HST 
would allow resolution of a spatial scale of '^1 kpc in 
the galaxy and measurement of a corresponding offset of 
the accreting SMBH from the galaxy core. The X-ray 
detection of SDSSJ0927-f-2943 holds promise for a deep 
Chandra study of SMBH and host galaxy. If the recoiling 
SMBH is a radio emitter, high- resolution radio observa- 
tions would provide even more precise measurement of 
its location. 

In summary, SDSSJ0927-I-2943 is the best candidate 
to date for a recoiling supermassive black hole. Its 
spectrum shows the characteristic signature of two sepa- 
rate emission-line systems, kinematically offset by ^2650 
km s~^ - broad lines from gas bound to the recoiling hole, 
and narrow lines from gas left behind. Further study 
of this source and detection of similar ones will provide 
important information on BH recoil, relevant timescales, 
and the frequency of BHs removed from their host galax- 
ies, and therefore on simulations of galaxy formation and 
evolution, and on the question how many BHs grew early 
by merging. 



We thank D. Merritt and G. Hasinger for illuminating 
discussions. We have made use of the SDSS, GALEX 
and ROSAT database. 
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